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Abstract 

The aim of this paper is to investigate discrete approximations of the 
exponential functional J °° exp(_B(t) — vt)dt of Brownian motion (which 
plays an important role in Asian options of financial mathematics) by the 
help of simple, symmetric random walks. In some applications the dis- 
crete model could be even more natural than the continuous one. The 
properties of the discrete exponential functional are rather different from 
the continuous one: typically its distribution is singular w.r.t. Lebesgue 
measure, all of its positive integer moments are finite and they charac- 
terize the distribution. On the other hand, using suitable random walk 
approximations to Brownian motion, the resulting discrete exponential 
functionals converge a.s. to the exponential functional of Brownian mo- 
tion, hence their limit distribution is the same as in the continuous case, 
namely, the one of the reciprocal of a gamma random variable, so ab- 
solutely continuous w.r.t. Lebesgue measure. This way we give a new, 
elementary proof for an earlier result by Dufresne and Yor as well. 



1 Introduction 

The geometric Brownian motion (originally introduced by the economist P. 
Samuelson in 1965) plays a fundamental role in the Black-Scholes theory of 
option pricing, modeling the price process of a stock. It can be explicitly given 
in terms of Brownian motion (BM) B as 

S{t) = S a exp (<rB(t) + (fj,- a 2 /2) t) , t > 0. 

In the case of Asian options one is interested in the average price process 

1 I* 

A(t) = - S(u)du, t > 0. 
t Jo 
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The following interesting result is true for the distribution of a closely related, 
widely investigated exponential functional of BM: 



2 = / exp(B(t) - vt) dt = — (v > 
Jo Z 2v 



0). (1) 



Here Zi v is a gamma distributed random variable with index 2v and param- 
eter 1, while = denotes equality in distribution. This result was proved by 



Dufresne (1990)] using discrete approximations with gamma distributed ran- 
dom variables and also by Yor (1992)] , using rather ingenious stochastic analy- 



sis tools. For more background information see [Yor (2001)] and [Csorgo, M. (1999)] 



As a consequence, the pth integer moment of I is finite iff p < 2v and 

E(2 P) = 2P r(2 ^- p) . (2) 

On the other hand, all negative integer moments, also given by @, are finite 
and they characterize the distribution of T. 

The situation is much nicer when BM with negative drift is replaced in the 
model by the negative of a subordinator (at,t > 0), that is, by the negative of 
a non-decreasing process with stationary and independent increments, starting 
from the origin. Then, as was shown by |Carmona et al. (1997)] , all positive 
integer moments of J — J Q °° exp(— a t ) dt are finite: 

nJP) = *(!)*!. *<p) > $(A) = -i l0 S E ( ex P(-^)), 0) 

and in this case the positive integer moments characterize the distribution of 
J. 

To achieve a similar favorable situation in the BM case, at least in an ap- 
proximate sense, it is a natural idea to use a simple, symmetric random walk 
(RW) as an approximation, with a large enough negative drift. Besides, in some 
applications a discrete model could be more natural than a continuous one. It 
seems important that, as we shall see below, the discrete case is rather different 
from the continuous case in many respects. 

So let (Xj)°^ 1 be an i.i.d. sequence with P(Xi = ±1) = \ and So = 0, 
Sk = Ylj=i {k> 1). Introduce the following approximation of I: 

oo 

Y = e MSk -kv) = l + £ 1 + £i£ 2 + • ■ • , & = exp(X J - i/), (4) 

k=0 

where v > 0. In this paper we investigate the properties of Y, which will 
be called the discrete exponential functional of the given RW, or shortly, the 
discrete exponential functional. 

In Section 2 below it turns out that the distribution of Y is singular w.r.t. 
Lebesgue measure if v > 1 . Then in Section 3 we find a formula, similar to (]2J , 
for all positive integer moments of Y when v > 1, and, because Y is bounded 
then, these moments really characterize its distribution. Finally, in Section 4 
we use a nested sequence of RWs to obtain a.s. converging approximations of 
I, and this way an elementary proof of result |T]) of Dufresne and Yor as well. 
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2 The distribution of the discrete exponential 
functional 



Let us start with a natural generalization: (£,j)°^ 1 be i.i.d., £j > 0. Consider 
first the finite polynomial 

Y n = 1 + + . .. + &...&, (5) 

= i + a(i + 6 + 66 + --- + 6---Cn) («>i), 

Yo = 1. This implies the following equality in distribution: 

r„ = i + (6) 

where £ = £i, and £ is independent of Y n -\. Since 7„/y = l+d +?i^2 H 

a.s., we get the basic self- similarity of Y in distribution: 

^ = 1 + (7) 

where £ is independent of Y. We remark that infinite polynomials similar to 
Y were studied by Vervaat (1979)] and many others. There some of the ideas 



discussed below have already appeared. 

A standard application of the strong law of large numbers gives a condition 
for having an a.s. finite limit Y here, see Theorem 1 in [Szekely, G. (1975)| . 
Namely, when E(|log£j|) < oo, one has Y n /• Y < oo a.s. if and only if 
E(log&)<0. 

In the special case when Y is defined as in but S n is the partial sum of 
an arbitrary i.i.d. sequence (Xj)°°^ 1 with zero expectation, Y < oo a.s. iff the 
drift added is negative: v > 0. Hence this condition is always assumed in our 
basic example (simple, symmetric RW). 

Next we want to show that self-similarity ([7} implies a simple functional 
equation for the distribution function F(y) — P(Y < y), y £ K. For a modest 
generalization of our basic case, let us introduce some notations. In ([5]) let £j 
take the positive values 71 < • • • < 7jv, and let pi = P(£ = 7,-). (In our basic 
case N = 2, 71 = e~ x ~ v , 72 = e 1_iy , pi = P2 = 5.) Consider the following 
similarity transformations: Ti(x) = 7^2; + 1 (1 < i < AT). When 

AT 



E(log£) = log 7l < 



holds, by © we have P(Y < y) = P(l + £Y < y) = J2i=iPi p (! + < 
y\£, = 7i) = Sj=i P» P(l + — u)- Thus one obtains the following functional 
equation for the distribution function of Y: 



N 



F(y) = J2PiF(T t r 1 (y)). (8) 
i=i 

An important special case is when 7^ < 1 (in the basic case: v > 1). 
Then by ([5]), Y is a bounded random variable. Moreover, each mapping Ti is a 
contraction, having a unique fixpoint = (1 — 7;)~ 1 , < i/i < •• • < yjv < 00. 



3 



Since each 7j is an increasing function, Ti(yj) < Tj(yfc) if j < k. Also, Tj(yi) < 
Tk(yi) if J < fc- Then it follows that each Ti maps the fundamental interval 
I = [3/1,2/ ;v] into itself. Clearly, I contains the range of Y too. 

In the case 7^ < 1 it is useful to rephrase the given problem in the language 
of fractal theory, see e.g. |Falconer (1990)1 . Let us introduce the symbolic space 
£ = {i = (ii , 12, . ■ ■ ) : ij = 1, . . . , -/V}, endowed with the countable power of the 
discrete measure (pi, ■ ■ ■ ,Pn), denoted by P. By ([5]), 

Y n = 1 + 7il (1 + 7<a (. • ■ (1 + 7i J)) = (T n o T J2 o • • • o T <B )(1), (9) 

with probability Pi t Pi 2 ■ ■ -Pi n - Thus the canonical projection II : £ — ► I, IL(i) = 
lim fc ^ 00 (T ll o • • • o T ik )(l) = limfc^.oo(l + 7i , + • • • + 7il . . . 7i J maps E onto 
the range of Y. The attractor A of the iterated function scheme of similarity 
transformations (Ti, . . . , Tn) is defined as 

N 

A=f| |J (T il o...oT i J(7), A=(jT i (A). (10) 

fe>0 l<h,...,i k <N i=l 

Then A is a non-empty, compact, self-similar set. In (|10p the fundamental 
interval I = [yi,yjv] can be replaced by any interval J which is mapped into 
itself by each T i; e.g. by J = [0,yjv] 3 1. Thus range(F) C A, cf. ©. 
The converse is also true, since for any y G A and for any e > there is an 
i G S and a large enough fc such that y G (X^ o ••• o Ti k )(J) and the length 
ICZfc o • • • o T ifc )(J)| = 7il ■ ■• 7 *JJ| < e. Hence, by ©, y G range(F), that is, 
A =rangc(y). Also, the distribution of Y on the real line, which will be denoted 
by Py, is simply P o II -1 . 

We are going to use the notations 

Jfi...i fc = [s/u...t„i,yii...i fc jv] = ( T »i "- oT *J( / )> ( n ) 

Vi!...i k l = ■ • • T ik )(yi) (I = l,...,N) as well, where ^ = 1, ... N and 
yi = (1 — 7 ;)~ 1 . The length of such an interval is [_Z^ 1 i fc = 7l • ■ ■ 7 i fc |/|, where 

|i1 = vn - yi- 

Returning to the distribution of Y in the basic case, consider first when the 
intervals h = Ti(I) = [3/11,3/12] = [2/1,2/12] and I 2 = T 2 (I) = [2/21,2/22] = [2/21,2/2] 
do not overlap, where 3/12 = 1 + 7 i(l — 7 2)~ 1 , 2/21 = 1 + 72(1 — 7i) ■ Thus 
there is no overlap iff yi 2 < 2/21, i-e., v > log(e + e~ r ) s» 1.127. Since F(y) = 
if y < 2/1 an d F(y) = 1 if y > 3/2, in this non-overlapping case © simplifies as 

f iF(Tf 1 (3/)) if ye [3/1,3/12), 

^(2/) = < 1 ifyG [yi2,y 2 i), (12) 

I i + H^GO) if2/G [y 2 i,y 2 ). 

By the similarities given by Ti and T 2 , applied to (TT2"j) . one obtains that 
-F has constant value | over the interval [ym, 3/121) an d constant value | on 
[y2i2,y22i)- Continuing this way by induction one gets that F has constant 
dyadic values over such plateau intervals: 

k 

F(y)-2- fc - 1 +^(z J -l)2-''', y G [y n ... 4fc i2,y 4l .., ifc 2i), i* = l,2. (13) 
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The sum of the lengths of these plateaus is \I\ (1 — (71 + 72)) (1 + (71 + 72) 

+(71 + 72) 2 H ), so add up to |/|. Hence the attractor A (the range of Y), 

i.e., the set of points of increase of F, has zero Lebesgue measure. So it is a 
Cantor-type set: an uncountable, perfect set of Lebesgue measure zero. 

The distribution function F is clearly a continuous singular function. For, 
if j/o G A and e > is given, take k so that 2~ fe < e. By the construction of A, 
there exists an interval Ii 1 ..,i k 3 Ho- Let the left endpoint of the left neighbor 
plateau of Ii t ...i k be rj\ (or —00), and the right endpoint of the right neighbor 
plateau be 772 (or 00). If 5 = min(?/o —T)i)V2~ Vo) > 0, then for any y such that 
\y - y \ < 5 one has \F(y) - F(y )\ < 2~ k < e by (gSD . 

It is not difficult to see, cf. [Grincevicius (1974)] , that in general, any solution 
of ([7|) has either absolutely continuous or continuous singular distribution. 

We mention that standard results of fractal theory, see Theorem 9.3 in 
Falconer (1990)] , imply that the Hausdorff dimension s of A equals the (unique) 



solution of the equation 7J 5 + 7J = 1. Solving this equation for i>, we get 
v = s _1 log(e s + e~ s ). Hence the fractal dimension s is a strictly decreasing 
function of v > log(e + e -1 ), tending to 1 as v — > log(e + e _1 ) and converging 
to as v — > 00. Also, the Hausdorff measure of A is H S (A) = |J| S , where the 
Hausdorff dimension s is the one defined above. It means that 



H S (A) 



^(l-e(2coshs)- 1/s ) 1 - (l-e- 1 (2coshs)- 1/s ) ^ 



Thus "H S (A) ->■ e 2 - e~ 2 as v ->• log(e + e _1 ) and H S (A) -> as v -)• 00. 

Next we are going to show that the distribution of y is singular w.r.t. 
Lebesgue measure even in the overlapping case if v > 1. Again, we con- 
sider the slight generalization introduced above. The proof below is based on 
Simon, K. et al. (2001)] and on personal communication with K. Simon. 



Theorem 1. Let £ take the values 7, {i = 1, . . . , N), < 71 < ■ • • < jn < L 
and let pi = P(£ = 7.;). Take an i.i.d. sequence (£j)i=ii £j = £• Then the 
distribution of Y = 1 + £1 + £1^2 + ■ ■ ■ is singular w.r.t. Lebesgue measure, if 



N N 



-Xp = E(logf) = ^Pilog7i < 5^ ^ log pi = -Hp. 

i=l i=l 

T/izs wi/Z &e called the entropy condition. Here xp is the Lyapunov exponent of 
the iterated function scheme {T\, . . . , T/y) corresponding to the Bernoulli mea- 
sure P. 

Proof. We are going to use the fractal theoretical approach and notations in- 
troduced above. 

We want to show that 

(DP Y )(x) = limsup P ^, {X,r J ) = 00 Pya.s, (14) 

where B(x, r) denotes the open ball (in the real line) with center at x and radius 
r and A is Lebesgue measure. The statement of the theorem easily follows from 
this. For, take the set E = {x e I : (DP Y )(x) = 00}. Then (Jl4j) implies that 
Py{E) = 1, while e.g. Theorem 8.6 in Rudin (1970)] shows that the symmetric 



derivative DPy exists and is finite A a.e., so X(E) = 0. 
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Introduce the notation ai (£) = : ii = j, 1 < I < k}. Thus 

oo N 

fe=ii=i 

By the SLLN, the set A 3 ■ = || € £ : k^a^ '(i) — > pj j has probability 1 for every 

j = l,...,N and so has A = f]f =1 A r Let C = {x e I : DT^a;) n A ^ 0}. 
Then Py(C) = 1. 

If x 6 C, there exists i <E A such that II(i) = x and & _1 aj^(i) — > p^ as 
k — y oo for all j — 1, . . . , TV. Fix such an i and x. Let r/c be the smallest radius 
such that B(x, r^) D 1^ i k , where i = (ii, . . . , . . . ) and 1^ i k is dehned by 

The following facts are clear: (a) x € A, moreover, a; S Ii 1 ..,i k i see (fit))) 
and (JTTJ) ; (b) § | Ji x ...i fc | < r k < clI^.^J, where c> 1 is arbitrary; (c) = 

M nf=i 7^ (d) Py(B( S ,^)) > Pr(4...,) - P(*i. • •.,**) = IljliP""® 
Using these facts it follows for any k > 1 that 

, , „ /rriV 

p y (%^» > f2d/n -i f n^iPj fc _ 

A(B(x,r fc )) - l 1 U \ U N fc-^o'O) 

\llj=i /j 

By our assumptions concerning a; and i, the ratio on the right hand side 
converges to (p^ 1 ■ • -p^T) / (7^ • • • 7^™) as fc — > 00. The entropy condition of 
the theorem implies that this latter ratio is larger than 1. Hence (Tl4l) holds, 
and this completes the proof. □ 

Returning to our basic case, consider the entropy condition when 71 = 
er 1 -" ', 72 = e 1 -" and p 1 = p 2 = \. The condition holds iff v > log2 « 0.693, 
since this is equivalent to |(— 1 — 2/) + |(1 — 1/) < | log 5 + 5 log |. Combining 
this with the condition 72 < 1, this means that the distribution of Y is singular 
w.r.t. Lebesgue measure for any v > 1. 

Characterization of the distribution of Y when < v < 1 remains open. In 
that case one of the two similarity mappings, T2, is not a contraction anymore, 
and that situation requires more sophisticated tools than the ones above. 




3 The moments of the discrete exponential func- 
tional 

Let us consider first the general case: (Cj)^Li i.i.d., £j > 0, as at the beginning 
of the previous section. Now we turn our attention to the moments of Y = 
1 + £1 + £1^2 -I ■ If Y n is defined by ([5]), the equality in law ((BJ) implies 

E(5?) = E ((1 + CYn^r) = J2 (fc) MfcE(yti), (15) 

where p > integer and /ifc = E(£ fc ). As n — > 00, by monotone convergence we 
obtain 

E(Yn=p(fyv k V(Y k ). (16) 
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In (fTS"]) and (TTrj|) both sides are either finite positive, or +00. 

Theorem 2. Lei be an i.i.d. sequence, £j > 0, and Y = l + £i + 

• • • . For p > 1 reaZ, E(Y P ) < 00 if and only if fi p = E(£ p ) < 1. Then \i q < 1 
/or any 1 < g < p and E(Y P ) < (1 — p p ^ p )~ p as well. In this case if p > 1 is an 
integer, we also have the recursion formula 

E(n = T ^E©^ ) . (17) 

Proof. First, ([6]) and the simple inequality (l + x) p > 1 + x p (x > 0, p > 1, real) 
imply that E(F^) > 1 + ^ p E(y r f_ 1 ). Suppose that fj, p > 1. Since E(Y" P ) = 1, 
taking limit as n — >• 00, one gets that F,(Y P ) = 00. 

Conversely, suppose that fj, p < 1. Then by Holder's (or by Jensen's) inequal- 
ity, fi q < p3j v < 1 for any 1 < q < p as well. We want to show that E(Y" P ) is 
finite. Let us begin by observing that E ((Yj - F J -_ 1 ) P ) = E ((^ • • • ^) p ) = fjp p 
(j > 1). Hence by the triangle inequality and Yq = 1 we get that 

n n 1 

e^) 1 ^ < i + E E - ^--i) p ) 1/p = E < — 17?' 

3=1 3=0 1 Mp 

for any n > 1 when /i p < 1. Taking limit as n — >■ 00, it follows that ~E(Y P ) < 
(1 - p}J p )- p < 00. 

Finally, the recursion formula (|17p directly follows from (|16p when /i p < 1, 
p > 1 integer. □ 

For integer p > 1 it follows from (fTTj) by induction that E(F P ) is a rational 
function of the moments fj,i, . . . , p p : 

where the coefficients of the numerator are universal constants, independent of 
the distribution of £j . 

These universal coefficients a!f^ „■ in the numerator of (IT81) make a sym- 
metrical, Pascal's triangle-like table if they are ordered according to multiindices 
{jp-i, ■ ■ ■ ,ji) as binary numbers, see the rows p — 1, ... ,5: 

1 

1 1 

12 2 1 

1 3 5 3 3 5 3 1 

14969 16 11 4 4 11 16 96941 

Rather interestingly, based on the above table, one may conjecture that each co- 
efficient 3 - % is equal to the number of permutations 7r of the set {1, . . . ,p} 
which have descent ir(i) > 7r(i + 1) exactly where ji = 1, 1 < i < p — 1. Another 
interesting conjecture is that the following recursion holds: 

*tX = E «£Li-x. (n,..., lp )E{o,ir, P >i, 

Cji,...,ip-i)eS(ii,...,i») 
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where S(i\, . . . ,i p ) is the set of all distinct binary sequences obtained from 
(ix, ...,i p ) by deleting exactly one digit; — 1. For example, a^ 10 = 11 = 
a i 10 + a oio + a oi l ■ This recursion would imply that the above table contains only 
positive integers and has the symmetries af i +1 \ — a\ p + ^ \ = 1 _ i — 

b+i) 

Ul l-ip,...,l— ii ■ 

There is a nice analogy between the moments of the exponential functional 
of a subordinator and the moments of Y, compare ([3]) and (TT8|) . First, the sum 
of the coefficients in the numerator of (TT5|) is pi, as can be seen by induction. 
For, if one explicitly writes down E(Y" P ), based on the recursion (fTT|). taking a 
common denominator, the numerator of each earlier term except the last one 
is multiplied by factors 1 — In the sum of the coefficients of the numerator 
it means a multiplication by zero. On the other hand, in the last term one 
multiplies the numerator of E(Y" P_1 ) by pfi p -i, which results the sum p\ of the 
coefficients by the induction. 

Second, there is a relationship between the denominators of © and (fT8"]) 
as well. In the special case when Y is defined as in (j?]), but S n is the partial 
sum of an arbitrary i.i.d. sequence (Xj)°^ 1 with zero expectation, $(A) = 
— n~ x logE (exp(A(S„ — vn))) = — logE(£ A ), so $(fc) = — log^fc, corresponding 
to the factors in the denominator of . The factors 1 — in the denominator 
of l|18p are tangents to these. 

Finally, let us consider the moments of Y in our basic case. Then jik = 
E(£ fc ) = exp(— kv) cosh(fc). Since cosh(/c) < e k for any k > 0, it follows that 
\ik < 1 for any k > 1 when v > 1, therefore all positive integer moments of Y 
are finite in this case by Theorem [2j In particular, in Section 2 we saw that 
Y is a bounded random variable if v > 1, hence the positive integer moments 
characterize its distribution. On the other hand, when < v < 1, only finitely 
many moments of Y are finite. For, /j,f. > 1 if < v < fc _1 logcosh(fc) /• I 
as fc — > oo. For example, even \i\ > 1 (and consequently all E(F P ) = oo) if 
< v < logcosh(l) w 0.43378. 



4 Approximation of the exponential functional 
of BM 

In this final section we are going to show that taking a suitable nested sequence 
of RWs the resulting sequence of discrete exponential functionals Q converges 
almost surely to the corresponding exponential functional T of BM. Based on 
this, using convergence of moments, we will give an elementary proof of theorem 
Jl]) of Dufresne and Yor. 

The underlying RW construction of BM was first introduce d by [Knight (1962)] , 
and simplified and somewhat improved by Revesz (1990)] and [Szabados (1996)] 



This construction starting from an independent sequence of RWs (S m (k),k > 
0)m=oi constructs a dependent sequence of RWs (S m (k), k > 0)^ =0 by "twist- 
ing" so that the shrunken and linearly interpolated sequence (B m (t) = 2~ m 
S m (t2 2m ),t > 0)^ =0 a.s. co nverges to BM (B (t),t > 0), uniformly on bounded 
intervals, see Theorem 3 in |Szabados (1996)] . 

We need one more result about this approximation here. This is stated in a 
somewhat sharper form than in the cited reference, but can be read easily from 
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the proof there. Namely, see Lemma 4 in Szabados (1996) , for almost every ui 
there exists an mo(w) such that for any to > mo(w) and for any K > e, one has 

sup sup \B m+j (t)-B m (t)\< (log K)im2-%. (19) 

j>l 0<t<K 

Lemma 1. Let B m (t) = 2~ m S m (t2 2m ), t > 0, m > 0, be a sequence of shrunken 
simple symmetric RWs that a.s. converges to BM (B(t),t > 0), uniformly on 
bounded intervals. Then for any v > 0, as m — > oo, 

oo 

* m ^ 



k=0 



exp (2- m S m {k) - vk2~ 

X = / exp (B(t) — vt) dt < oo a. 
Jo 



Proof. The basic idea of the proof is that the sequence of functions f m (t,u>) = 
exp (B m (t) — vt), converges to f{t,uj) — exp (B(t) — vt) for t 6 [0, oo) as m — >• 
oo, for almost every lo. If one can find a function g(t,uj) € i 1 [0, oo), that 
dominates each f m for m > mo(w), then their integrals on [0, oo) also converge 
to the integral of /, and then we are practically done. 

First, by (|19|) . for a.e. w there exists an tuq — mo(w) so that for any K > e, 

sup sup \B m (t)-B mo {t)\<Ki (log K)i <K?logK, (20) 

m>m 0<t<K 

where we supposed that mo was chosen large enough so that mo2 _m °/ 2 < 1. 
Second, by the law of iterated logarithms, 

,. B mo (t) S mo (u) 
hmsup r = hmsup- 1 



t^oo (2tloglogi)2 m^oo (2u log log it) 2 

where u = t2 2m ". Hence for a.e. uj, there is a Kq = Kq(v,lu), such that for any 
t > K , 

B mo (t) < 2(floglogt)4 < 2ti logt, (21) 

where Kq is chosen so large that 3t? logt < vt/2 for any t > Kq. 

Since a.s. any path of B mo is continuous, it is bounded on the interval 

[0, Kq\. Then by (|20p , we have an upper bound uniform in to: for any m > too 

and t e [0,K ], B m (t) < M(uj). On the other hand, when t > K , by (j2~Tj) . 

B mo (t) < 2t? logt and so by (|2"0"]) . B m (t) < 3t? logt, for any to > mo- 
Summarizing, the function 

e MM if 0<t <K (y,u), 
e-^ 2 i£t>K (v,u), 

is an integrable function on [0,oo), dominating exp(_B m (i) — vt) for each to > 
toq(o;). This implies that 



/•oo />oo 

lim / exp (B m (t) — vt) dt = exp (B 
■"^°°Jo Jo 



(t) — vt) dt < oo a.s. 

Finally, compare / °° exp (B m (t) - vt) dt to Y m = 2- 2m J2T =0 ex P (B m (k2- 2m ) 
-vk2~ 2m ) that appears in the statement of the lemma. Applying the uniform 
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domination of exp (B m (t) — vt) by the function g shown above, both the tail of 
the integral on the interval [Kq, oo) and the tail of the sum for k > \Ko2 2m ~\ 
is smaller than J K exp(— vt/2) dt, thus their difference is uniformly arbitrarily 
small for any to > too if Kq is large enough. On the interval [0, -Kq] the dif- 
ference of the integral and the sum (which is a Riemann sum of a continuous 
function) tends to zero uniformly as m — > oo, since on each subinterval of length 
2" 2m , the difference of B m (t) and B m (k2- 2m ) is at most 2~ m . This completes 
the proof of the lemma. □ 

Next we want to apply the results of the previous sections to Y m . To do this 
we introduce the following notations. For m > and n > 1 let 

n 

Y m , n = 2- 2m ^exp(2-™5 m (fc)-^2- 2m ) 

fc=0 

_ 2 2m (1 -)- £ ml -|- £ ml £ m2 + (ml • • • Cmn) (22) 

and Y mfi j= 2- 2m , where £ mj = cxp(2-™X m (j) - v2-^). Here X m (j) = 
S m (j) ~ S rn(j - 1) (j = 1, 2, . . .) is an i.i.d. sequence, P(X m (j) = ±1) = |. 

Then F m ,„ F m as n — ► oo, F m < oo a.s. iff ^ > 0, and 1^ satisfies the 
following self-similarity in distribution: 

— 2 -)- ^ m Y 7n Or 2 = 1 + (m^ ^jth (23) 

where £ m and Y m are independent, £ m = £ m j. Using the notations of Section 2, 
now 71 = exp(-2-" l -j/2- 2m ), 72 = exp(2" m -i/2- 2m ), pi = p 2 = \. \iv > 2 m , 
72 < 1 holds, so the similarity transformations T\ and T 2 are contractions, 
mapping the interval / = [(1 — 7 i) _1 , (1 — 72) _1 ] into itself. By Theorem [TJ 
the distribution of Y m is singular w.r.t. Lebesgue measure if v > 2 2m log2 
(to > 1). Moreover, there is no overlap in the ranges of T\ and T 2 iff v > 
2 2m log(2 cosh(2~ m )). As to — ► oo this means asymptotically that v > | + 
2 2m log2 + o(l). 

For to > and fc integer let 

fi mk = nO = exp(-isk2- 2m ) cosh(fc2-" 1 ). 

Since Theorem [2] is applicable to 2 2m Y m , one obtains that E(y,^) < oo if and 
only if fi mp < 1 and then the following recursion is valid for p > 1 integer: 

E(l£) = f^2- 2m ^ fc ) Mmfe E(^). (24) 

Now using cosh(x) < e x (x > 0), it follows that fi m k < exp(fc2 _m (l — i>2~ m )). 
So p, m k < 1 for any k > 1 if v > 2 m . If < v < 2 m , only finitely many positive 
moments are finite, since n m k > 1 when < v < 2 2m k~ 1 logcosh(/c2~ m ) — > 2 rn 
as k — > oo. 

More importantly, 

fi mk < exp (k2- 2m Q - iX\ (k > 1), (25) 
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since cosh(ai) < exp(ir 2 /2) when x > (compare the Taylor series). Thus 
Hmk < 1 for any to > if v > ~. This condition is sharp as to — » oo. For, apply 
e x = 1 + x + o(x) and cosh(a;) = 1 + x 2 /2 + o{x 2 ) (as x — > 0) to the definition 
of [i mk . Then 

Hmk = l + k2- 2m (^-A+ o(2- 2m ), (26) 
for any fixed k as to — » oo. 

Lemma 2. If p is a positive integer such that § < ^, tten 

lim E(Y£) = ——J j- < oo. (27) 

Proof. By ()25[) . for any positive integer p such that | < ^ we have /i mp < 1. 
Since Theorem [2] is valid for 2 2m Y m , the recursion formula (|24| holds, and by 
induction one gets E(l^) as a rational function of the moments /i m i, . . . , /i m p, 
similarly to formula (fT5|) . The argument below formula (IT81 also applies here 
too, showing that the sum of the coefficients in the numerator of this rational 
function is p\2~ 2mp . The extra factor comes from the difference that Y m is 
multiplied by 2 2m here, compare equations (|17l) and ([24|. Since each n m k — > 1 
as to — y oo, it follows that 2 2mp times the numerator tends to pi. 

By dUD, we get that 1 - Mmfc = k2~ 2m {v - §) + o(2~ 2m ) if fc is fixed and 
to — > oo. So 2 2mp times the denominator of the rational function tends to 
pi Y\£—i(v — f ) as to — )• oo. This and the limit of the numerator together imply 
the statement of the lemma. □ 



Our next objective is to give an asymptotic formula, similar to (|27[). for the 
negative moments of Y m as to — > oo. 



Lemma 3. for o/Z integer p > I, we have 

lim E(Y-*) - lim E(y- 1 ) JJ ( v + ^ ) , (28) 



' k 



fc=l 



w/iere linim^oo E(F m 1 ) < oo. 



Proof. We want to show (j28|) by establishing a recursion. Introduce the no- 
tations z, m ^ k = E(Y~ k ) and fJ, m ,-k = E(£~ fe ) for k > 1 integer. By (|22J|, 
< Y^ 1 < 2 2m , hence all negative moments z mj & of F m are finite. 

The self-similarity equation (|2"3"|) implies that £,rnYm = Y m — 2~ 2m and so 



t-ly-l A f 



i - 2- 2 ™y ro - 1 ' 



where £ m and Y m are independent. Taking feth moment (k > 1 integer) on both 
sides and applying the Taylor series 



x k x (n — \ 



{\-x) 1 
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valid for any |x| < 1, one obtains 

n=k ^ ' 

with the notations introduced above. This implies that 

(/i m> _fe - l)z m ,k - k2~ 2m z nhk+ i = a(m, k), (29) 

where 

a(m,k) = £ ( n k Z%- 2m ^~ k h min >0. 

n=k+2 ^ ' 

Next we want to give an upper bound for a(m,k), which goes to zero fast 
enough as m — > oo. Since £ m > 71 = exp(— 2 _m — ^2 _2m ), by (f2"2"j) it follows 
that 

y-l < 2 2m A _ 2 2m^ _ ^ < 2 2m (2"" + ^2~ 2 " 1 ) < 2" l+1 , 

if to > log(f)/log(2), where we used that 1 — e~ x < x, for any real x. This 
implies that z m , r+i = E(r m ^) < E(y- r ) (sup^- 1 )) 3 ' < z m ^ m+1 ^ for 
r i j > 0. Substituting this into the definition of a(m, A;), one gets that 

a(m,k) < z„ hk+ i £ f n -^\ 2 - 2 ^- k h^+^ n - k -^ 



n=k+2 



= Zm , k+1 2 k ^-v-™-i jr ( n k _% [1 - m)n 

n=k+2 ^ ' 

= z^k+12-" 1 - 1 ((1 - 2 1 - m y k - 1 - k2 1 ~ m ) 

< z mik+1 2- m -Hk(k + l)2- 2m = z m , k +i2k(k + 1)2~ 3 "\ 

if to is large enough, depending on k. 

Let us substitute this estimate of a(m, k) into (|29l) and express the following 
ratio: 

Zrn,k+\ ^m, — k 1 

z m , k = fc2- 2 ™(l + 0(2-™))' 
Apply the asymptotics (|2l)l) here with — k: 

z m .k+\ = ^+|+o(l) 
z m ,h ~ 1 + 0(2-™) ' 

as to — > 00. This implies the equality lim m _ i . 00 z m< k+i/ z m ,k — v + f ; an d thus 
for any positive integer p, 

lim S?=fT^+-V OO) 

It remains to show that E(l r ^ 1 ) has a finite limit as to — > oo. Writing Y,^ 2 = 
y m l^" 3 and applying the Cauchy-Schwarz inequality, one obtains (E(F m 2 )) < 
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E(^)E(F m 6 ), or E(F m 2 ) < E(^)E(r m 6 )/E(y m 2 ). Suppose first that v > 1 
and take limit here on the right hand side as rn — > oo, applying (f2~Tf and ([3"U)) . 
It follows that sup m>1 E(F^ 2 ) < oo. As E(F^ 2 ) is an increasing function of 
v by its definition, hence the same is true for any v £ (0, 1] as well. Since by 
Lemma [TJ Y m —> X a.s., where each Y m and also I take values in (0, oo) a.s., 
it follows that Y^ 1 — > I^ 1 a.s. Then by the L? uniform boundedness of Y^ 1 
(m > 0) shown above, ~E(Y~ 1 ) — > E(/ _1 ) < oo follows as well. This ends the 
proof of the lemma. □ 

Finally, it turns out that Y^ 1 converges to I -1 in any LP . This makes it 
possible to recover the result ([T]) of |Dufresne (1990)] and [Yor (1992)] . 



Theorem 3. Let B m (t) = 2- ,n S m {t2 2m ), t > 0, m > 0, be a sequence of 
shrunken simple symmetric RWs that a.s. converges to BM (B(t),t > 0), uni- 
formly on bounded intervals. Take 

oo 

Y m = 2- 2m Y exp (B m (k2- 2m ) - vk^- 2m ) = 2~ 2m (1 + £ m i + + • ■ • ) 



k=0 



and 



f °° 

X = exp (B(t) - vt) dt 
Jo 



when v > 0. Then the following statements hold true: 

(a) Y~ converges to X^ 1 in LP for any p > 1 real and lim m _ ) . 00 ~E(Y rn p ) = 
E(X-f ) < 00; 

(b) X = 2/Z2U, where is a gamma distributed random variable with index 
2v and parameter 1; 

(c) Y m converges toX in LP for any integer p such that 1 < p < 2v (supposing 
v > ^) and then linim^oo E(l^) = E(I P ) < 00. The same is true for any real 
q, 1 < q < p. 

Proof. By Lemma U Y m — > X a.s., where each Y m and also X take values in 
(0, 00) a.s. Hence Y~ 1 — > X^ 1 a.s. By LemmaO lim m _ i . 00 E(Y~ k ) < 00 for any 
k > 1 integer, so (a) follows. 

Thus by (a) and Lemma [3J for any integer p > 1, 

p- 1 / . \ 

a p = E(2T p ) = c Y[ [v + -= ) = c2 1 - p (2v + 1) • • • (2v + p - 1) < 00, 



k=l 



where c = E(I 1 ). By a classical result, see Simon, B. (1998)] , a Stieltjes 



moment problem is determinate, that is the moments uniquely determine a 
probability distribution on [0, 00), if there exist constants C > and R > such 
that a p < CR p (2p)l for any p > 1 integer. In the present case a p < c2~ p (p+ 1)! 
when v < 1 and a p < c(v/2) p ~ 1 (p+ 1)! when v > 1, so the moment problem for 
I -1 is determinate and it also follows that I -1 has a finite moment generating 
function in a neighborhood of the origin. 

Also, using the moments of the gamma distribution we get 

b p = m- p Z p 2u ) = 2- ^^^ = 2^(2^(2*. + 1) • ■ • (21/ +p - 1), 
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for any p > 1, and Zivll has a finite moment generating function in a neighbor- 
hood of the origin as well. Writing down the two moment generating functions 
by the help of the moments a p and b p , respectively, it follows that 

E(exp(ttX -1 )) = -E(exp(uZ 2 „/2)) 

in a neighborhood of the origin. Substituting u = 0, one obtains c = E(Z -1 ) = v 
and this proves (b). 

Finally, again, Y m — > X a.s. by Lemma [1] If p is an integer such that 
1 < p < 2v, by Lemma [2J using the moments of the gamma distribution, and 
by (b), we have lim™^ E(Y r P) = E(2PZ^) = E(27). This proves (c). □ 
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